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Abstract 

A somewhat unorthodox method for determining vertical crustal motion at a tide-gauge location is to difference the sea level time 
series with an equivalent time series determined from satellite altimetry. To the extent that both instruments measure an identical ocean 
signal, the difference will be dominated by vertical land motion at the gauge. We revisit this technique by analyzing sea level signals at 28 
tide gauges that are colocated with DORIS geodetic stations. Comparisons of altimeter-gauge vertical rates with DORIS rates yield a 
median difference of 1.8 mm yr~‘ and a weighted root-mean-square difference of 2.7 mm yr” 1 . The latter suggests that our uncertainty 
estimates, which are primarily based on an assumed AR(1) noise process in all time series, underestimates the true errors. Several sources 
of additional error are discussed, including possible scale errors in the terrestrial reference frame to which altimeter-gauge rates are 
mostly insensitive. One of our stations, Male, Maldives, which has been the subject of some uninformed arguments about sea-level rise, 
is found to have almost no vertical motion, and thus is vulnerable to rising sea levels. 

Published by Elsevier Ltd. on behalf of COSPAR. 
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1. Introduction 

The precise measurement of vertical crustal motion is 
one of the most important, yet highly challenging, tasks 
of modern geodesy, with applications to a wide range of 
geophysical problems. The signals are tiny — of order 1 cm 
per year or smaller — even for “large” motions such as those 
associated with Pleistocene ice unloading. The present 
capability of space geodetic methods to measure vertical 
motion is discussed by (among others) Altamimi et al. 
(2007), Collileux et al. (2007) and Riilke et al. (2008). 

Determination of vertical land motion at tide gauges is 
fundamental to the measurement of sea level, be it by tide 
gauges or by satellite altimetry (Douglas et al., 2001; 
Woodworth, 2006). For tide-gauge measurements of abso- 
lute sea level, knowledge of vertical motion is obviously 
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critical. For altimetry the requirement arises through cali- 
bration/validation, which is often performed via a global 
network of tide gauges; insufficiently known vertical 
motion at the gauges is currently the limiting factor in 
accurate altimeter validations (Mitchum, 2000). For these 
reasons the international community has organized large 
geodetic campaigns that are beginning to provide the nec- 
essary data (Schone et al., 2009; Woppelmann et al., 2009). 

Another approach to recover land motion at tide gauges 
is to turn the calibration problem around and difference the 
common sea level time series seen by both the altimeter and 
the gauges. To the extent the sea level signals are truly com- 
mon, the difference is a measure of vertical land motion at 
the gauge. The approach has been explored by Cazenave 
et al. (1999), Nerem and Mitchum (2002), Garcta et al. 
(2007), and others. We revisit the problem here because 
as the time series lengthen we expect some level of 
improved precision in all estimates, and we now have in 
hand over sixteen years of high-quality altimetry from 
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the Topex/Poseidon (T/P), Jason-1, and Jason-2 satellites. 
Moreover, the data processing issues are sufficiently com- 
plex to warrant several approaches to the problem, and 
our approach, while fundamentally similar to previous 
work, is different in some important details. Finally, as part 
of this special journal issue devoted to the DORIS geodetic 
system, we take the opportunity to compare recent 
DORIS-derived vertical rates with our own altimeter-gauge 
estimates. 

2. Data processing 

In briefest terms our approach is to form near-global 
gridded arrays of altimetric sea levels every 9.92 days (the 
repeat period of the satellites) and to extract from these 
grids sea-level time series near selected tide-gauge loca- 
tions. The tide gauge data are low-pass filtered to corre- 
spond to the signals in the gridded altimetry. Differencing 
any two such time series gives an estimate of the vertical 
land motion at the gauge location. The following subsec- 
tions give further details for each data type. 

2.1. Altimeter data 

We use data from the satellites Topex/Poseidon (Fu 
et al., 1994), Jason-1 (Menard et al., 2003), and Jason-2. 
All data have been collected along the primary T/P ground 
track from September 1 992 through March 2009. This time 
series of altimetry comprises 608 repeat cycles, although 
some repeats have only partial or no data because of oper- 
ational issues. In our time series the transition dates 
between the three satellites are 4 May 2002 and 12 July 
2008. Our analysis begins with official project-released ver- 
sions of the data as follows: T/P from MDGR-B, Jason-1 
from GDR-C, and Jason-2 from standards consistent with 
the Jason-1 GDR-C. We have then updated a number of 
altimeter correction fields, the most important of which 
we note below. 

A great deal of effort has been expended to ensure the 
data are as consistent as possible between the different sat- 
ellites and that both the altimeter ranges and the myriad 
required corrections are well calibrated and validated. 
Beckley et al. (2004) describe in some detail the work 
involved in seamlessly merging T/P and Jason-1 data, 
and similar work has recently been done for Jason- 1 and 
Jason-2 (Beckley et al., 2010). The cross-satellite transitions 
have been markedly helped by the special tandem-mission 
calibration campaigns during which the ‘old' and ‘new’ sat- 
ellites fly in close proximity for several months. On the 
basis of these near-simultaneous ocean elevation measure- 
ments we find a 99 mm bias between T/P and Jason- 1 
and a 175 mm bias between T/P and Jason-2. (Note that 
these bias estimates are to some extent functions of range 
corrections, orbits, and other variables; see, for example, 
Dorandeu et al. (2004).) 

Within-mission calibrations are also critical. For exam- 
ple, the onboard microwave radiometers, necessary to cor- 


rect altimetric ranges for wet tropospheric delay, have been 
found prone to drift (Keihm et al., 2000; Ruf, 2000) and are 
a continuing source of attention (Scharroo et al., 2004; 
Brown et al., 2004; Zlotnicki and Desai, 2004). We use a 
recalibrated Topex microwave radiometer replacement 
product to MGDR-B (Brown et al., 2009) and a Jason- 1 
microwave radiometer replacement product for GDR-C 
(Brown et al., 2007). See the cited references for details of 
the calibration techniques. 

We use a completely new set of orbital ephemerides 
computed for all three satellites by Lemoine et al. (submit- 
ted for publication), based on DORIS and laser tracking of 
the satellites with ground-based tracking stations situated 
according to the International Terrestrial Reference Frame 
2005 (Altamimi et al., 2007). The frame to which our esti- 
mated vertical crustal motions refer is a somewhat subtle 
issue. The altimetric satellites, of course, orbit the center 
of mass of the whole Earth system and not the center of 
frame formed by the tracking stations. The orbit determi- 
nation process partially compensates for the difference, 
mostly through empirically derived accelerations adjusted 
to fit the tracking data. The amount of compensation, how- 
ever, is not always certain; it depends in part on modeling 
issues involving the geocenter, a topic of current research in 
orbit determination. In any event, through the tie to the 
tracking stations the altimetric sea-surface heights will tend 
to be approximately located within the ITRF2005. The ref- 
erence frame used for the tracking stations does clearly 
affect the final sea-surface heights (Morel and Willis, 
2002; Morel and Willis, 2005; Beckley et al., 2007), espe- 
cially in the Z (or northward) direction, but also in orthog- 
onal directions as well. To the extent that the altimetric sea- 
surface heights are located in the ITRF2005, we consider 
our estimated vertical station motions are also approxi- 
mately in that frame. 

There is, however, one important exception. Altimetric 
sea levels are mostly immune to scale errors in the reference 
frame. This was pointed out by Morel and Willis (2005) by 
experimenting with scale-dependent adjustments to the 
geodetic stations tracking T/P, and it has been confirmed 
in unpublished orbit-determination calculations by our 
own group. Changing the reference-frame scale may inflate 
(or deflate) the satellite tracking residuals, but it otherwise 
leaves the semimajor axis — and hence mean sea level — 
unaffected, the axis being fixed primarily by the satellite 
orbital period. Thus, our altimeter-gauge vertical land 
motions are also insensitive to reference-frame scale, in 
contrast to more direct determinations by geodetic mea- 
surements. We return to this point below in the context 
of comparisons with DORIS. 

Because the issue of residual drift in the altimetry is of 
critical importance to inferred rates of land motion, any 
additional possible monitoring or estimation of drift is 
worthwhile. We have therefore subjected our final altimeter 
time series to the tide-gauge validation procedures 
described by Mitchum (2000). Prof. Mitchum kindly 
computed the data shown in Fig. 1 , which shows mean 
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Fig. 1. The results of the altimeter validation procedure computed 
courtesy of Prof. Gary Mitchum on our final time series of altimetric 
sea-surface heights. Each dot gives the mean difference (mm) between 
altimetric and tide-gauge sea levels evaluated at a large array of tide 
gauges. See Mitchum (2000) for further details, including the procedures 
used to account for land motion at tide gauges. TOP diagram shows the 
outputs before we have corrected for inter-satellite bias. BOTTOM 
diagram shows after correction. Colors denote satellite according to: T/P 
(cyan), Jason-1 (dark blue), Jason-2 (purple). 


differences between altimetric and tide-gauge sea levels 
computed for every 10-d repeat cycle; the top diagram 
shows the mean differences before correcting for the 99- 
mm and 175-mm inter-satellite biases, the bottom diagram 
shows the mean differences after correcting. A linear fit to 
the bottom data gives a trend of only —0.02 mm yr -1 . The 
Mitchum procedure itself is sensitive to reference frame 
uncertainties and site instabilities, as Mitchum has stressed, 
but to the extent that the procedure is a valid test for sys- 
tem drift, our data show no drift and moreover show no 
evidence of any residual bias between missions. 

It is important to note that tide gauge comparisons, such 
as Fig. 1 , are used only to test and validate the altimetric 
heights, not to calibrate them. The altimetry and gauge 
data remain independent in our analysis. 

To form 10-day gridded altimeter sea-surface heights we 
have found it useful to apply a correction for atmospheric 
loading, since it removes substantial power at periods 
below 10 days and thereby reduces noise and improves 
the consistency of neighboring tracks during the gridding 
step. For this we used an inverted barometer correction, 
based on daily mean pressures from NCEP Reanalysis-2 
(Kanamitsu et al„ 2002) weather products. Daily means 


are used rather than the more standard 3-b or 6-h products 
because (a) the ocean is manifestly not isostatic at sub-daily 
frequencies and (b) S? atmospheric tide signals in the 
higher-frequency pressure data interfere with the tidal cor- 
rection of altimetry, the models for which already properly 
account for S 2 radiational tides (Ray, 1999). 

2.2. Differencing altimetry and tide gauge data 

All tide gauge data used in this paper are from the 
archives of the University of Hawaii Sea Level Center. 
For the main analysis part of the paper, we selected tide 
gauge stations situated within 150 km of a DORIS geodetic 
station, although we also use a few additional gauges to 
illustrate certain points. We also require the overlap 
between the altimeter and gauge time series to be at least 
five years. In our experience shorter segments yield unreli- 
able vertical rates. There are 28 stations satisfying these cri- 
teria; their positions are shown in Fig. 2. 

We use daily mean sea levels, which are computed at the 
Hawaii center from hourly data by means of a 119-point 
convolution filter which has a response of 95, 50, and 5% 
at periods of 124, 60, and 40 hours. This essentially elimi- 
nates all tidal power in the diurnal and higher frequency 
bands. 

Since we estimate vertical land motion by subtracting 
the sea levels observed by altimetry and tide gauges, it is 
essential that the processing of both data types be as con- 
sistent as possible. For example, we must correct the tide- 
gauge data for the tides of long period (weekly through 
18.6 yrs), plus the pole tide, with models consistent with 
those models used for the altimetry. The details of the mod- 
els need not concern us, nor even their accuracies, since 
common errors will cancel; it is the consistency of models 
that must be ensured. Perhaps the most important is allow- 
ance for the 18.6-yr node tide, for which we use a self-con- 
sistent equilibrium model (Agnew and Farrell, 1978), which 
can be roughly approximated by a degree-2 zonal spherical 
harmonic with amplitudes approximately 8 mm along the 
equator and double that at the poles. Neglecting to apply 
this correction to the tide gauges could, if we are unlucky, 
induce in a 9-year time series a false trend of about 2 mm 
yr -1 , or even double that in high latitudes. 

Likewise, since an inverted barometer correction was 
used for altimetry, we applied the identical correction 
(based on mean daily atmospheric pressures) to the daily 
tide-gauge data. 

Finally, to attempt to match the temporal characteristics 
of the gridded altimeter data we applied a low-pass filter to 
all tide-gauge data. The filter has a (broad) cutoff at a per- 
iod around 25 days. 

2.3. DORIS vertical rates 

A number of groups have processed DORIS data and 
computed station positions and velocities which could be 
used to compare against our altimeter-gauge vertical rates. 
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Fig. 2. Locations of tide gauges and DORIS stations used in this study. Unfilled circles are tide gauges referred to in the text but for various reasons (e.g., 
no nearby DORIS data) are not included in Tabic 1. 


We prefer, however, to use time series of station coordi- 
nates, since these can sometimes shed more light on station 
motion. From a time series of station heights we can com- 
pute velocity trends as needed. 

In this paper we have employed the recent series solution 
of P. Willis and colleagues at the Institut Geographique 
National (IGN), denoted ‘ign09wd01 .stcd’ which consists 
of weekly estimates of station coordinates (Willis et al., 
2010) available from the International DORIS Service 
(Tavernier et al., 2006). The ign09wd01 solution is a deriv- 
ative of the solution submitted by the IGN for ITRF2008, 
and includes updates to the orbit processing for radiation 
pressure modelling and drag parameterization that 
improve the station position, geocenter and earth-orienta- 
tion determination. The IGN series solutions are also 
essentially equivalent to those earlier described by 
Gobinddass et al. (2009). DORIS solutions have been pro- 
duced by other groups as well; three solutions that contrib- 
uted to the ITRF2005 have been independently assessed by 
Altamimi et al. (2006). Of the seven DORIS solutions 
recently submitted to ITRF2008, the Willis et al. (2010) 
solution is at present the only one available as station-coor- 
dinate time series. 

From these station coordinate time series we compute 
vertical velocities, with each weekly point weighted by its 
formal error. Problematic data flagged by Willis et al. 
(2009, Table 1) are not used. One delicate point is how to 
handle multiple occupations at a given site. In such cases 
it is possible to merge time series by using local geodetic 
ties, but we have not attempted to do so. Instead we simply 
compute vertical velocities for each individual occupation 
and compute a weighted mean of these for our final station 
velocity estimate. This is emphasized in our time series 
plots below (e.g.. Figs. 9-12) where we display each sepa- 
rate DORIS occupation period with zero mean eleva- 
tion — i.e,, the time series will purposely appear disjointed. 


2.4. Error bars 

All of the altimeter-gauge time series that we have tested 
fail the widely used Durbin-Watson test for serial correla- 
tion (Draper and Smith, 1998), usually with very high sig- 
nificance levels (i.e., more strictly, they fail the null test of 
no serial correlation). In light of this, throughout this paper 
whenever we estimate a linear trend from a time series, we 
account for serial correlation by assuming an order- 1 auto- 
regressive AR( 1 ) noise process with the correlation coeffi- 
cient computed directly from residuals. Such an error 
model is often adopted in geophysical time series analysis, 
since red noise processes are so common (Bloomfield and 
Nychka, 1992, e.g.; von Storch and Zwiers, 1999). A vari- 
ety of other noise models for DORIS station-coordinate 
time series have been extensively explored by Williams 
and Willis (2006), including a model equivalent to AR(1); 
in fact, the autoregressive model, combined with a variable 
white noise model (dependent on weekly formal errors in 
DORIS estimates), was preferred in at least one of their 
tests, although in the end the authors recommended a mod- 
ified flicker noise model. Some of their error models (espe- 
cially the exponential decaying white noise model) are 
unlikely to apply directly to altimeter and tide-gauge time 
series, and it is outside the scope of the present paper to 
investigate in detail the most applicable noise model. Our 
approach here is to use the AR(1) model as a useful general 
approximation for handling serial correlation. Lee and 
Lund (2004) give a brief review of the problem and give 
useful formulae for confidence intervals which we adopt 
here. 

For the altimeter-gauge vertical rates we attempt to 
account for one further systematic error, which is the 
uncertainty in the overall altimeter drift rate as deduced 
from a Mitchum-type calibration/validation study (Nerem 
and Mitchum, 2002, cf. ). This latter uncertainty was 




1514 


R.D. Ray et al. i Advances in Space Research 45 (2010) 1510-1522 


thought by Mitchum (2000) to be roughly 0.4 mm yr” 1 , 
stemming primarily from uncertainties in the land motion 
at gauges. This is generally smaller (often much smaller) 
than the autoregressive-based standard errors on our esti- 
mated vertical rates, but it implies that our altimeter-gauge 
rates have error bars never less than this adopted value of 
0.4 mm yr” 1 . 

There is a simiar systematic error in our DORIS-based 
vertical rates which we do not explicitly account for. 
Because the DORIS weekly positions have already been 
transformed into the ITRF2005 via 7-parameter transfor- 
mations, the quoted DORIS uncertainties are in some sense 
internal measures of those transformations without regard 
to the underlying errors in the ITRF itself (P. Willis, pers. 
commun., 2010). For vertical rates the largest errors 
induced by reference-frame errors arise from uncertainties 
in the global Z (northward) drift of the frame and in the 
scale of the frame. Since many of our tide-gauge stations 
are in low to middle latitudes, the scale may be the more 
important. Present-day uncertainty in the scale is thought 
to be of order 0.1 ppb yr” 1 (Z. Altamimi, pers. commun., 
2008), which translates to about 0.6 mm yr” 1 in the verti- 
cal. This is often larger than our quoted standard errors 
in DORIS rates, so our quoted errors should be considered 
somewhat too optimistic. Note that any error in the ITRF 
scale would affect all DORIS-based vertical rates by pre- 
cisely the same amount. 

All quoted uncertainties below represent l-a standard 
errors. 

3. Revisiting Cazenave et al. (1999) 

Before turning to a more comprehensive compilation 
and comparison with DORIS measurements, we revisit a 
few interesting examples of altimeter-gauge vertical motion 


estimates which were originally discussed by Cazenave 
et al. (1999). Those authors called attention to several loca- 
tions displaying anomalously large vertical rates. The most 
extraordinary was at Socorro Island, off the western coast 
of Mexico, where time series of altimeter-gauge differences 
and weekly DORIS coordinate solutions both showed a 
very large ~30 cm subsidence over a span of about three 
years. This apparently reflected deformation following a 
volcanic eruption in 1993. Unfortunately, we can add little 
to their discussion because the tide gauge measurements at 
Socorro discontinued in Sept. 1997. However, DORIS 
measurements have continued, and they form the basis 
for a recent geophysical analysis of Socorro Island by Bri- 
olc et al. (2009). 

Of the other locations of anomalous vertical motion 
noted by Cazenave et al. (1999), the three largest rates were 
at Saipan, Nuku Hiva, and Rabaul. None of these three 
sites had or has DORIS measurements. In Fig. 3 we show 
an updated time series for Saipan. It is clear that the large 
vertical motion of 20.0 ± 2.1 mm yr” 1 noted by Cazenave 
et al. abruptly stopped just after their paper was published. 
Since then there has been a slow apparent subsidence, 
which mostly cancels the earlier rise; a linear fit to the 
entire time series gives a trend of 0.6 ± 0.6 mm yr” 1 . 
Cazenave et al. stated that the cause of the observed trend 
at Saipan was “unclear from a geological point of view,” 
and we concur, especially in light of the observed change 
in trend. There is no record of any instrument changes 
around 1997-1998. 

At Nuku Hiva the large negative rate noted by Cazenave 
et al. (1999) appears to have continued (Fig. 4), although 
the tide gauge time series was interrupted over the interval 
1998-2004. At this point we are unsure about the relative 
datums before and after the interruption, so we have split 
the figure. Nonetheless, estimated linear trends for both 



Fig. 3. (Bottom) Allimetric and tide-gauge sea levels at Saipan. Atmospheric loading signals have been removed from both time series by use of an 
inverted barometer correction based on daily mean atmospheric pressures. (Top) Differences between altimctric and tide-gauge sea level time series, each 
point representing a 10-d estimate. To the extent that the altimeter and the tide gauge observes identical sea-surface elevations, then the difference 
(altimeter minus gauge) is a proxy for vertical position of the tide gauge. Cazenave et al. (1999) called attention to the high vertical rate at Saipan during 
the 1992-1997 time period, but that rate abruptly stops thereafter. Fitted linear trend (magenta line) is for the entire time series. 
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Fig. 4. Altimetric and tide-gauge sea levels at Nuku Hiva (bottom) and their differences (top), a proxy for vertical position of the tide gauge. Tide gauge 
measurements were interrupted over the interval 1998-2004. 


sections of the time series are comparable, but if we assume 
continuity we find a trend for the entire series significantly 
lower: -7.1 ±2.7 mm yr 1 . Moreover, the early section 
might be better explained as nearly constant (zero trend) 
but with one discontinuity sometime when the time series 
was interrupted near beginning 1995; the later section after 
2004 is more complicated. 

At Rabaul we again have no further tide gauge data 
than those employed by Cazenave et al. (1999). However, 
in this case our Fig. 5 is not so much an update to their fig- 
ure as it is a correction in light of better documentation. As 
Cazenave et al. noted, a series of volcanic events in late 
1994 caused significant subsidence at the Rabaul tide 
gauge. This information was reported by the Rabaul 
Volcanological Observatory to the Hawaii Sea Level Cen- 




RED <■ Ode gauge; BLUE = altimetry 


1994 1996 1998 

Fig. 5. Altimetric and tide-gauge sea levels at Rabaul (bottom) and their 
differences (top), a proxy for vertical station position. The large discon- 
tinuity in early October 1994, representing subsidence of over 15 cm, 
followed a series of nearby volcanic eruptions. Large, roughly linear, 
subsidence continues for several more years until the tide-gauge data stop. 


ter and the tide gauge data were subsequently adjusted by 
— 13 cm in October 1994 and by —14 cm in November 1994 
and thereafter (Patrick Caldwell, pers. commun., 2009). 
Thus, the Rabaul inferred land motion in Fig. 6 of Cazen- 
ave et al. (1999) is much smoother than it should have 
been. If we “uncorrect” the data (our Fig. 5) we see the 
large jump (~15 cm) reported by the Rabaul personnel. 
In addition, a large, roughly linear, subsidence continues 
for several more years until the tide-gauge data stop. 

In the spirit of Cazenave et al. (1999), we here report a 
recent unusual implied land motion (or possible tide-gauge 
problem) at Colombo, Sri Lanka — see Fig. 6. The altime- 
ter-gauge differences display a large jump sometime in the 
spring of 2006. To our knowledge there were no significant 
earthquakes in that region during that time, so the time 
series may well indicate some instrumental problem. The 
DORIS time series at Colombo was discontinued in 



Fig, 6. Altimetric and tide-gauge sea levels at Colombo, Sri Lanka 
(bottom) and their differences (top). 
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October 2004, owing to siting difficulties, so DORIS data 
can shed no light on this issue. 

In the time series data for our 28 primary stations we 
have noticed a few cases of similar changes in datum level, 
although none so dramatic as at Colombo. To the extent 
that these can be reliably recognized we try to account 
for them in any trend estimates. Our experience suggests 
(see also discussion below) that further efforts to unravel 
these possible datum problems, as well as to improve 
tide-gauge documentation, would be a worthy endeavor. 

4. Estimated vertical rates and comparisons with DORIS 

The estimated vertical rates at our 28 selected tide-gauge 
stations and at the corresponding DORIS stations are tab- 
ulated in Table 1. The altimeter-gauge and DORIS rates 
are more directly compared in Fig. 7. From the table and 
the diagram it is readily apparent that the altimeter-gauge 


error bars tend to be larger than the DORIS error bars. 
The median error bar is 1.30 mm yr” 1 for the altimeter- 
gauge data, but 0.54 mm yr” 1 for the DORIS data. The 
spread in the altimeter-gauge vertical rates is also greater, 
perhaps reflecting the larger error bars. By happenstance 
most DORIS rates are small, and only three stations have 
rates near or exceeding 5 mm yr” 1 . 

The dotted line in Fig. 7 represents a linear lit to all 28 
points. Although the error bars in the DORIS estimates are 
mostly small, they are not zero, so we have used an ‘errors- 
in-variables’ type regression that allows for errors in both 
variables (Madansky, 1959; Lybanon, 1984). The slope is 
0.92 ± 0.11 and the intercept is — 1.12 ± 0.26 mm yr -1 
(an ordinary least squares regression that ignores errors 
in the DORIS data yields a biased slope of only 0.60). 
The weighted coefficient of correlation is 0.56, where we 
have weighted each point by the reciprocal sum of squares 
of the individual x and y error bars. The significance of this 


Table 1 

Estimated vertical crustal motions from altimeter-gauge and DORIS (ign09wd01) data. 


UH 

id 

GLOSS 

id 

Tide Gauge 

Gauge 

Location 


C.I. 

(%) 

RMS 

(cm) 

H All -Gauge 

(mm/y) 

DORIS 

Sites 

Dist. 

(km) 

No. 

wks 

H DORIS 

(mm/y) 

AH 

(mm/y) 

002 

113 

Betio, Tarawa 

1.4° N 

172.9° E 

96 

2.2 

0.0 ±0.5 

betb 

1.1 

127 

—2.4 ±1.7 

2.4 

003 

169 

Baltra, Galapagos 

0.4° S 

90.3° W 

80 

2.3 

0.5 ±0.6 

serb 

130.2 

194 

2.3 ±1.0 

-1.8 

015 

140 

Papeete, Polynes. 

17.5° S 

149.6° W 

97 

2.1 

—0.7 ± 0.5 

papb paqb 
path 

7.1 

661 

0.1 ±0.5 

-0.8 

016 

138 

Rikitea, Polynes. 

23.1° S 

134.9° W 

94 

3.8 

-1.2± 1.0 

rikb 

1.3 

112 

2.7 ±2.2 

-3.8 

019 

123 

Noumea, N. Caledonia 

22.3° S 

166.4° E 

65 

3.8 

4.8 ±1.8 

noua nowb 

3.7 

548 

1 .2 ± 0.6 

3.8 

022 

137 

Easter Is. 

27.2° S 

109.4° W 

95 

6.6 

-2.2 ±3.0 

easa easb 

6.4 

769 

2.3 ±0.3 

-4.5 

053 

149 

Guam 

13.4° N 

144.6° E 

85 

3.2 

0.9 ±0.7 

guab 

30.9 

320 

2.7 ±0.8 

-1.8 

058 

_ 

Kauai, Hawaii 

22.0° N 

159.4° W 

100 

3.8 

-0.2 ±1.0 

koka kolb 

36.9 

849 

1.1 ±0.4 

-1.2 

079 

128 

Chatham Is. 

43.9° S 

176.6° W 

41 

2.6 

-2.9 ± 1.1 

chab 

111.0 

441 

0.7 ±0.3 

-3.6 

081 

175 

Valparaiso, Chile 

33.0° S 

71.6° W 

97 

7.2 

-4.2 ± 1.7 

sanb 

91.0 

413 

-2.1 ±0.3 

-2.1 

094 

- 

Matarani, Peru 

17.0° S 

72.1° W 

91 

9.0 

— 13.0 ±2.7 

area arfb 

89.0 

583 

-6.1 ±0.7 

-6.9 

108 

028 

Male, Maldives 

4.2° N 

73.5° E 

99 

2.1 

— 1.0 ± 0.5 

malb 

1.2 

183 

-0.8 ±1.0 

-0.2 

121 

273 

Pt. LaRuc, Seychelles 

4.7° S 

55.5° E 

99 

3.1 

-4.9 ± 0.5 

mahb 

0.8 

344 

-3.9 ±0.8 

-1.0 

161 

_ 

Jakarta 

6.1° S 

106.8° E 

57 

6.6 

1.3 ±2.5 

cibb cicb 

41.4 

738 

— 1 .6 ± 0.4 

2.9 

179 

024 

Amsterdam Is. 

38.7° S 

77.5° E 

70 

6.4 

-8.9 ±2.6 

amsa amtb 
amub 

101.4 

497 

4.9 ±0.6 

-13.9 

180 

023 

Kerguelen 

49.4° S 

70.2° E 

97 

3.6 

4.3 ±0.8 

kera kerb 
kesb ketb 

2.7 

751 

3.6 ±0.5 

0.6 

211 

245 

Pt. Delgada, Azores 

37.7° N 

25.7° W 

69 

3.6 

-l.0± 1.7 

pdlb pdmb 

1.6 

468 

2.0 ±0.5 

-3.0 

223 

253 

Dakar, Senegal 

14.7° N 

17.4° W 

52 

3.5 

0.7 ± 1.3 

daka 

6.3 

382 

1.4 ±0.6 

-0.7 

235 

329 

Palmeira, C. Verde 

16.8° N 

23.0° W 

41 

3.0 

-2.3 ±2.8 

salb 

5.5 

263 

3.2 ± 0.6 

-5.5 

276 

223 

St. John's 

47.6° N 

52.7° W 

95 

8.0 

3.0 ± 2.1 

stjb 

4.3 

458 

0.4 ±0.3 

2.6 

288 

229 

Reykjavik 

64.2° N 

21.9° W 

42 

3.2 

-2.4 ± 1.3 

rcya reyb 
rezb 

3.0 

797 

-0.8 ±0.3 

-1.7 

291 

263 

Ascension Is. 

7.9° S 

14.4° W 

39 

2.0 

0.1 ±0.9 

asdb 

9.3 

492 

-2.2 ±0.5 

2.3 

292 

264 

St. Helena 

15.9° S 

5.7° W 

64 

1.7 

0.3 ±0.6 

hela helb 
hemb 

6.0 

620 

1.7 ±0.5 

-1.5 

370 

073 

Manila 

14.6° N 

121.0° E 

58 

5.8 

-7.5 ±3.2 

mana manb 

9.6 

694 

1.8 ±0.5 

-9.3 

554 

159 

La Jolla 

32.9° N 

117.3° W 

93 

4.8 

-2.1 ± 1.4 

monb moob 

78.3 

187 

2.9 ± 1.9 

-5.0 

569 

- 

San Diego 

32.7° N 

117.2° W 

100 

4.7 

-2.3 ± 1.2 

monb moob 

73.0 

187 

2.9 ± 1.9 

-5.2 

600 

181 

Ushuaia, Argen. 

54.8° S 

68.3° W 

42 

6.1 

3.9 ±4.2 

rioa riob 
ripb riqb 

118.8 

816 

5.2 ± 0.5 

-1.8 

755 

218 

Virginia Key, FL 

25.7° N 

80.2° W 

81 

5.8 

-3.3 ±2.2 

miab 

0.2 

228 

2.0 ± 0.7 

-5.2 


C.I. — Completeness Index, denoting percent of gauge data available in 1993-2009 (after any edits). 
RMS = root-mean-square difference between altimetry and gauge data. 

Dist. = distance between DORIS station and tide gauge. 

No. wks — number of weekly DORIS solutions available at this site. 

AH — difference in altimeter-gauge and DORIS vertical rates. 
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Fig. 7. Scatter diagram comparing estimated vertical land motion rates 
from our altimeter-gauge analyses with those inferred from DORIS 
solution ‘ign09wd01.stcd’. The dotted line is an ‘errors-in-variables' 
regression fit to the data, with slope 0.92 ±0.11 and intercept -1.12 ± 
0.26 mm yr" 1 . Amsterdam Island is seen to be a case of particularly poor 
agreement between the two kinds of vertical-rate estimates; see Fig. 10. 

correlation exceeds 99% according to results of a permuta- 
tion test, a standard resampling technique for testing statis- 
tical significance (Good, 2009). 

For all 28 stations the RMS difference between altime- 
ter-gauge and DORIS vertical rates, weighted according 
to the individual station standard error estimates, is 2.7 
mm yr" 1 . The median difference is —1.8 mm yr" 1 . The 
reduced x 2 statistic is 3.6, large enough to suggest that 
our vertical-rate standard errors are underestimated, a 
point anticipated in our discussion in Section 2.4. In fact, 
of the 28 differences between altimeter-gauge and DORIS 
vertical rates, 19 exceed the range of the associated quoted 
errors, far more than should be expected. If we arbitrarily 
inflate all error bars by 50%, then 10 exceed the range, 
which is closer to the expected number for 1 -a standard 
errors. 

A close examination of Table 1 reveals that several of 
the stations showing poor agreement between altimeter- 
gauge and DORIS rates should not be considered truly 
colocated. For example, the DORIS station at Arequipa, 
Peru, is located 89 km from the tide gauge at Matarani, 
but more importantly it is in the Andes at elevation 
2450 m. In such a tectonically active region we cannot 
expect any consistency in vertical motion. Indeed, in the 
case of Arequipa we have the additional difficulty of a large 
earthquake occurring on 23 June 2001, which produced a 
vertical offset of several cm and at least a temporary change 
in velocity (Melbourne et al., 2002; Willis et al., 2009). (For 
this reason we had already removed the ‘areb’ occupation, 
corresponding to the period 2002.0-2003.6, from our tabu- 
lation; the pre-earthquake ‘area’ and post-earthquake 
‘arfb’ occupations do yield consistent vertical rates.) 


Similar points can be made for the distance between 
Monument Peak and the tide gauges at San Diego and 
La Jolla. 

If we limit the tide gauge and DORIS station separation 
to 6 km we find many of the discrepant stations in Fig. 7 
drop out. Fig. 8 shows the revised scatter plot, now with 
14 points. The slope of the fitted line is 1.11 ± 0.16 and 
the intercept is —0.78 ± 0.39. The weighted coefficient of 
correlation is now increased to 0.78. For these 14 stations 
the RMS difference between altimeter-gauge and DORIS 
estimates drops to 2.0 mm yr" 1 and the median difference 
is —1.2 mm yr" 1 . 

Although these kinds of statistical summaries are useful, 
it is enlightening to examine the individual time series. 
Space precludes a presentation of all 28, so we instead 
focus on four examples. Of these, two correspond to signif- 
icantly non-zero vertical rates, one being a ‘good’ example 
(it lies approximately along the unit-slope line of Fig. 7) 
and the other a ‘bad’ example (far from the line). These 
are shown in Figs. 9 and 10, respectively. 

The Amsterdam Island example (Fig. 10) shows medio- 
cre agreement between the altimetric and tide-gauge sea 
levels, with a correlation coefficient of only 0.50 and an 
RMS difference of 6.4 cm. There is also a strong suggestion 
of a change in slope near the data gap in 1999. Clearly if 
the altimeter and tide gauge are failing to observe a com- 
mon sea-level signal, then any inference about vertical land 
motion is suspect. That appears to be the case here. 

Nonetheless, the altimeter-gauge differences in Fig. 10 
do suggest a large negative trend, so the disagreement with 
the DORIS-based trend of +4.9 ± 0.6 mm yr" 1 is all the 
more striking. Note, however, the DORIS estimate is dom- 
inated mostly by the second occupation (amtb); the rate 
from more recent data (amub) is much smaller. An early 



Fig. 8. As in Fig. 7, except using only those points where the DORIS 
ground station is within 6 km of the tide gauge. The dotted line is an 
‘errors-in-variables’ regression fit to the data, with slope 1.1 1 ± 0.16 and 
intercept -0.78 ± 0.39. 
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1994 1996 1998 2000 2002 2004 2006 2008 2010 


Fig. 9. (Top) Altimetric and tide-gauge sea levels at Point La Rue in the Seychelles. The correlation between the two curves is 0.92. (Middle) Ten-day 
differences between the two time series shown in the top panel, a proxy for vertical position of the tide-gauge station. (Bottom) Time series of weekly 
vertical coordinates from DORIS solution ign09wd01 .sted. The precision of these estimates increases with time as more DORIS satellites become 
available. A least-squares fit to the DORIS data gives a linear trend of —3.92 ± 0.77 mm yr -1 . 



1994 1996 1998 2000 2002 2004 2006 2008 2010 


Fig. 10. As in Fig. 9, but for Amsterdam Island in the South Indian Ocean. Different colors in the bottom panel denote different DORIS occupations; each 
occupation is separately demeaned, so apparent discontinuities between different occupations may not be real. The DORIS data yield a linear trend of 
+4.93 ±0.55 mm yr' 1 , based primarily on occupation amtb (magenta color). The DORIS and altimeter-gauge vertical-rate estimates are in poor 
agreement, possibly owing to the poor agreement between altimetric and tide-gauge sea levels. 

DORIS solution by Soudarin et al. (1999) found a vertical fewer DORIS satellites than those now available. More- 
rate of —12.5 mm yr -1 . Although this is much closer to our over, Willis et al. (2009) point out that the Amsterdam 

altimeter-gauge rate, it is surely obsolete, being based on DORIS station has been affected by several antenna tilt 
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instabilities. The time periods of these antenna problems 
are not precisely known, but it does imply that some peri- 
ods reflect antenna motion rather than land motion. This 
problem was not discovered until after the Soudarin et al. 
(1999) study, so antenna motion may have contributed to 
their large negative rate. Overall, however, we attribute 
the poor agreement between DORIS and altimeter-gauge 
rates at Amsterdam Island to the apparent inability of 
our altimeter signal to reproduce the gauge signal (or 
vice-versa). Mitchum (2000) also points to “poor ocean sig- 
nal agreement” at this location, and for this reason he rou- 
tinely deletes this tide gauge from his altimeter validation 
tests. 

Kerguelen Island (Fig. 1 1 ) is another problematic site. 
Mitchum (2000) again points to a “possible tide gauge level 
shift,” From our altimeter-gauge difference time series 
(middle panel) we agree that around 1998 there is a fairly 
clear discontinuity in level, and also in trend. If we fit a 
trend to only the post- 1998 data (which we have done in 
Table 1) then we find good agreement with the DORIS- 
based rate of +3.6 ± 0.5 mm yr . 

The tide gauges at both Amsterdam Island and Kergue- 
len Island are in fact bottom-pressure recorders, which are 
prone to drift over time (e.g., Woodworth et al., 1996). The 
Kerguelen deployments are described in some detail by 
Testut et al. (2006), and those authors summarize various 
efforts made to determine the drift in the instruments. Over 
the period shown in our Fig. 1 1 , four successive pressure 
sensors have been installed in the same stilling well. The 
first was installed in 1993 and replaced in 1998, which is 


the point in time we detect a change in slope. Other replace- 
ments occurred in 2001 and 2004. The difficulties in recov- 
ering system offsets and drifts are emphasized by Testut 
et al. As our figure suggests, it may well be that satellite 
altimetry can help in these efforts. 

Finally in Fig. 12 we show data from Male in the Mal- 
dives Islands, a chain of low-lying coral islands southwest 
of India. As can be seen, the agreement between altimetry 
and tide gauge is very good; the RMS difference is 2.1 cm 
and the correlation coefficient is 0.89. The altimeter-gauge 
and DORIS vertical land rates are in good agreement, both 
near zero but negative, suggesting either no vertical move- 
ment or very slight subsidence. These results thus add fur- 
ther weight to the evidence, summarized by Woodworth 
(2005), that the Maldives are unlikely to escape the conse- 
quences of ongoing and future sea level rise. That there 
should be any question in the matter stems from some ill- 
informed arguments made mostly, but not completely, in 
the popular media (see criticisms by Kench et al., 2005). 
It seems clear that the low-lying Maldives are in a vulner- 
able position regarding sea level rise. 

5. Discussion and conclusions 

As Cazenave et al. (1999) and Nerem and Mitchum 
(2002) noted previously, the use of combined satellite 
altimeter and tide-gauge data offers a valid and encourag- 
ing approach to the problem of determining vertical 
motion at tide-gauge sites. While the more straightforward 
geodetic approaches based on DORIS and GPS network 




Fig. 11. As in Fig. 10, but for Kerguelen Island. The fitted trend (red line) to the altimeter-gauge differences is based on only the post-1998 data. The 
DORIS data yield a linear trend of +3.63 ± 0.51 mm yr -1 . Each DORIS occupation is separately demeaned when shown in the bottom panel. 
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Fig. 12. As in Fig. 9, but for Male in the Maldives Islands, southwest of India. A least-squares fit to the DORIS data gives a linear trend of —0.8 ± 1.0 mm 
yr -1 . Thus, both altimeter-gauge and DORIS estimates of vertical motion suggest that the Maldives Islands are stable over this time period, or even 
slightly subsiding. 


solutions are becoming more mature, the consistencies of 
those solutions are not always satisfactory. The altimeter- 
gauge approach may thus add some useful information 
to these analyses. 

As noted above, potentially useful information could 
stem from the general independence of the altimeter-gauge 
rates from scale errors in the terrestrial reference frame. 
Both the intercept estimates from Figs. 7 and 8 and the 
median differences between altimeter-gauge and DORIS 
vertical rates suggest a systematic global difference of 
roughly 1 ± 0.5 mm yr“ 1 , with DORIS measuring higher 
rates. A rather large error of order 0.2 ppb yr“ ! in the scale 
of ITRF2005 could explain this, but we are not inclined to 
push such a conclusion. Our large uncertainty, sometimes 
stemming from inconsistent altimetric and tide-gauge sea- 
level curves, the relatively large scatter in Figs. 7 and 8, 
and problems such as those highlighted for Amsterdam 
Island, all urge caution. Precisions in the altimeter-gauge 
vertical rates need to be improved before more definite 
statements can be made, but they clearly show potential. 

It is not surprising that the best altimeter-gauge rate esti- 
mates — best in terms of standard error of the estimate- as 
well as agreement with DORIS — are generally for those 
tide gauges on small open-ocean islands (Easter Island 
and Amsterdam Island are the exceptions to this). The 
island data show the best agreement between altimetric 
and tide-gauge sea level signals; see the RMS values listed 
in Table 1. Consistency in these two signals is, of course, 
desirable since the whole technique relies on the altimetry 
and the tide gauge seeing a common ocean signal. 


It is also not surprising that the largest RMS differences 
between altimetry and gauges occur for those stations 
along the coasts of large land masses. In contrast to the 
deep ocean, shelf and near-coastal ocean processes tend 
to display much shorter wavelengths, a fact that almost 
ensures the altimeter and gauge are observing somewhat 
different signals, especially since at some gauges the dis- 
tance to the nearest altimeter track can be over 200 km. 
Moreover, the altimetry can be impacted by a number of 
problems near the coast, including land contamination in 
the altimeter footprint and, more commonly, in the radi- 
ometer footprint and inaccurate environmental correc- 
tions, including tides. 

Fortunately, this suggests a path toward some progress, 
since the improvement of shallow-water satellite altimetry 
is an ongoing concern of several research groups (Vignud- 
elli et al., 2005, e.g.; Bouffard et al., 2008; Vignudelli et al., 
2010). Some of this work involves merging data from mul- 
tiple satellites to obtain the higher spatial resolutions 
needed in coastal zones. This is clearly desirable, but it will 
require more complex validation techniques to ensure the 
integrity of long-period signals, since some satellite altime- 
ters are designed for studying mesoscale phenomena, not 
subtle climatic signals (or vertical motion signals, as in 
the present case). In any event, this ongoing work aims 
to improve all satellite data through improved waveform 
processing and better corrections, so it should be immedi- 
ately applicable to studies such as ours. 

In addition to the study of crustal motion, the compar- 
ison of satellite altimeter and tide-gauge data has the 
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potential to reveal unknown problems in the latter. We 
noted several cases where the gauge data are either prob- 
lematic or insufficiently documented. See aso the discussion 
by Mitchum (2000) in this regard. 

Vertical crustal motion is one of the most challenging of 
geodetic measurements. Comparisons of new vertical rates 
from colocated GPS, DORIS, and other geodetic systems is 
one of the most important aspects of improving the terres- 
trial reference frame. More comprehensive and updated 
analyses of satellite altimetry and tide gauge data are war- 
ranted for the insight they can lend to these problems. 
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